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Abstract 


A zonal mean chemistry transport model (2-D CTM) coupled with a semi-spectral 
dynamical model is used to simulate the distributions of trace gases in the present day 
atmosphere. The zonal-mean and eddy equations for the velocity and the geopotential 
height are solved in the semi-spectral dynamical model. The residual mean circulation is 
derived from these dynamical variables and used to advect the chemical species in the 2- 
D CTM. Based on a linearized wave transport equation, the eddy diffusion coefficients 
for chemical tracers are expressed in terms of the amplitude, frequency and growth rate 
of dynamical waves; local chemical loss rates; and a time constant parameterizing small 
scale mixing. The contributions to eddy flux are from the time varying wave amplitude 
(transient eddy), chemical reactions (chemical eddy) and small scale mixing. In spite of 
the high truncation in the dynamical module (only three longest waves are resolved), the 
model has simulated many observed characteristics of stratospheric dynamics and 
distribution of chemical species including ozone. Compared with the values commonly 
used in 2-D CTMs, the eddy diffusion coefficients for chemical species calculated in this 
model are smaller, especially in the subtropics. It is also found that the chemical eddy 
diffusion has only a small effects in determining the distribution of most slow species, 
including ozone in the stratosphere. 
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1. Introduction 


Two-dimensional (zonal mean) chemistry transport models (CTMs) have played 
an important role in studying various physical and chemical processes which control the 
distributions of chemical constituents in the atmosphere. Comprehensive three- 
dimensional model with full chemistry is computationally too intensive to be used for 
long term (order of decades) simulation of the atmosphere. In addition, a zonal-mean 
approach provides an important alternative to interpreting the enormous amount of data 
from three-dimensional model and satellite observations. Therefore, the continuing 
improvements of two-dimensional models are indispensable. 

This paper documents a highly truncated semi-spectral dynamical model coupled 
with a two-dimensional photochemistry transport model of the stratosphere. With small 
increases in computing resources, the approach described here provides the following 
improvements from traditional two-dimensional models in the calculation of dynamical 
variables, the transport of trace gases, and the interactions between the dynamics and 
tracer distribution: 

• The primitive equations for the dynamical variables are separated into two sets, 
one for zonal-mean variables and one for eddies. Here, the word "eddy" refers to 
quantities that deviate from the zonal-mean. The eddy fluxes in the zonal mean 
momentum and thermodynamical equations represent transport of momentum and heat in 
the latitude-height plane due to zonally asymmetric motions and are calculated from 
resolved waves instead of using parameterization in terms of zonal-mean quantities. The 
three longest planetary waves are generated by imposing observed seasonally varying 
geopotential height at the lower boundary. Their propagation through the troposphere 
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into the stratosphere is calculated by the model. Kinnersley (1995) demonstrated that 
including wave-wave interactions is important in properly distributing energy among 
waves. Also included in the eddy equations is a linear damping term which represents 
the wave cascading to small unresolved scales. These small scale waves contribute to 
eddy mixing of chemical tracers when they are coupled with the planetary waves. 

• The zonal-mean concentrations of the trace gases are obtained by solving the 
zonal-mean transport equation. Transport in the latitude-height plane of a 2-D CTM 
consists of the residual circulation and the eddy flux that represents net transport due to 
zonally asymmetric motions. The residual circulation is calculated directly according to 
its definition from the Eulerian mean velocity and eddy flux of temperature. The 
contribution of eddy fluxes to the zonal mean distribution of chemical tracers in the 2-D 
CTM module is parameterized as a diffusion flux with the diffusion coefficients 
calculated explicitly from the dynamical waves. The approach of calculating eddy 
diffusion coefficients is based on the linearized eddy tracer transport equation similar to 
that used by Pyle and Rogers (1980), Strobel (1981), and Smith and Brasseur (1990). 
With the contribution from small scale mixing and chemical eddy added to the transient 
eddy, the total eddy transport is mainly diffusive (down-gradient). The values of eddy 
diffusion coefficients are smaller than those commonly used in 2-D CTMs especially in 
the tropics. 

• The model includes interactions between dynamics, radiation and chemistry. The 
chemical composition of the atmosphere can affect the dynamics by changing the 
heating/cooling rates through radiative transfer. Meanwhile, the dynamics can change 
the tracer distribution through changes in the residual mean circulation, zonal mean 
temperature and eddy diffusion coefficients. 
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The paper is organized as follows. Section 2 is a simple review of the eddy 
parameterization in 2-D CTMs. In Section 3 we give a brief description of the model, 
with emphasis on the dynamical module. The scheme of calculating eddy diffusion 
coefficients is formulated in Section 4. The model is used to simulate the present day 
atmosphere and the results for dynamical fields, the eddy diffusion coefficients and the 
distributions for trace gases, such as N 2 O, CH 4 and O 3 , are compared with observations 
and with previous studies in Section 5. Discussions of results are given in the last 
section. 


2. Parameterization of Eddy Fluxes 

The transformed Eulerian mean (TEM) circulation in the pressure coordinates or 
the zonal-mean circulation in the isentropic coordinates have been extensively used in 
two-dimensional chemistry transport models. If the zonal asymmetric motion can be 
described by small amplitude steady adiabatic waves, these two circulations become 
identical and represent the Lagrangian mean motions of air parcels (Tung, 1984). The 
total transport of trace gases is represented by an advection due to the circulations and an 
eddy flux which is mainly diffusive. The circulations are driven by diabatic heating, 
dissipation and eddy forcing. 

While it is commonly accepted that the eddy transport in the middle atmosphere is 
mainly due to organized planetary scale waves, rather than chaotic or turbulent mixing, 
there are still controversies about how to parameterize the eddy fluxes. A series of efforts 
to understand the mechanism of eddy transport and to derive the mixing coefficients 
began with Reed and German (1965). They parameterized the eddy fluxes in a two- 
dimensional model with Eulerian mean circulation as the product of a 2x2 diffusive 
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tensor and the gradient vector of zonal mean mixing ratio of tracers. However, due to an 
oversimplified assumption about the trajectories of air parcels, their results were later 
challenged. Matsuno (1980) showed that the trajectory of an air parcel driven by a steady 
planetary wave is elliptic when projected onto the meridional plane, and, therefore, there 
is no net transport for steady waves. Plumb (1979) and Pyle and Rogers (1980) reached 
the same conclusion using the linearized eddy tracer transport equation. Their results 
were the manifestation in Eulerian framework of the earlier study by Andrews and 
McIntyre (1978), which used the generalized Lagrangian mean formulation. 

During the past several years, using more powerful computers and new numerical 
techniques, the problem of transport by planetary waves has been revisited by following 
trajectories of a large number of air parcels in the stratosphere or using the ’contour 
advection technique' (Schoeberl and Bacmeister, 1993; Bowman, 1993; Waugh and 
Plumb, 1994; Chen, 1994). These studies have provided intuitive pictures of eddy 
mixing. However, due to the lack of small scale mixing, these methods are useful only 
for simulations of a short time period. It is difficult to define directly from their results a 
parameterization of eddy mixing suitable for a global 2-D CTM. 

The eddy flux for an inert tracer is given by the zonal mean of V • x ' , where v' 
is the eddy velocity in the meridional direction and x' * s the eddy tracer distribution. 

For steady waves, v' and %’ are orthogonal (90° phase lag between two waves) and the 
zonal-mean of their product, the eddy flux, is equal to zero. When the wave is not steady, 
the amplitude of %' changes with time. In this case, the phase relation between v' and 
X' changes due to the tendency term and/or chemical sink term in the wave equation and 
thus, generate net eddy flux. The eddy flux generated this way is related to the wave 
transience (transient eddy. Plumb, 1979) and chemical reactions (chemical eddy, Smith 
and Brasseur, 1990). Calculation of the transient eddy diffusion coefficient based on the 
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tracer wave equation have met with some conceptual and practical difficulties. It is easy 
to show that growing waves generate down-gradient eddy flux, and so are diffusive. On 
the other hand, decaying waves produce up-gradient eddy flux, which can initiate 
numerical instability in model calculations. As decaying and growing phases occur 
alternatively in the atmosphere, they would cancel each other. However, diagnostic 
calculations using observed wind suggest that the eddy transport is almost always 
diffusive (e. g. Plumb and Mahlman, 1987). 

A popular parameterization scheme uses the flux-gradient relation of potential 

v' * q ' 

vorticity (PV) to derive eddy diffusion coefficients ( K q = — ) either from 

dq/dy 

observed or from model-calculated (three-dimensional) wind fields (Tung, 1987; 

Hitchman and Brasseur, 1988; Newman et al., 1988; Yang et al., 1991). The derived 
eddy diffusion coefficient for PV, K yy , is then applied to all chemical species. As will 

be shown in Section 4, the eddy diffusion coefficient for any tracer (including PV which 

is a dynamical tracer) has a dynamical component (which is the same for all tracers), and 
a chemical component (which is species dependent). Therefore, the derived K yy are 

appropriate for other tracers only if the chemical eddy component of the eddy diffusion 
are small for PV and for chemical tracers. It is quite possible that using K yy for 

chemical species could overestimate the eddy transport for long-lived species by 
including the 'chemical eddy' of PV. There are other difficulties in practice. K q 

calculated from observed winds will contain negative values at some time and locations. 

These have to be ignored and substituted with positive values in a model simulation. 
Typical values of calculated K q are as large as 10 6 m 2 /sec. When these large values 

are used in two-dimensional models, they tend to underestimate the meridian gradients 
for long-lived gases in the stratosphere. 


6 



Garcia (1991) suggested a new method to parameterize eddy mixing through 
planetary wave breaking in the stratosphere. The time dependent equation for the eddy 
quasi-geostrophic potential vorticity ( q' ) is solved explicitly for a single planetary 
wave. The eddy equation contains dissipation terms representing damping by Rayleigh 
friction, Newtonian cooling and a separate linear term, -8q' , representing damping due 
to wave-breaking. The wave -breaking criterion is |<?y|/ ^ ^1. The value for the damping 

rate, 8 , is calculated at every time step for each location as follows. It is set to zero if 
|^|/5y <1* 1° region where >1, the values for 8 are calculated from the wave 

activity equation to stop further increase of the wave activity density. Garcia showed that 
the eddy flux for PV, v'q' , can be written in the form -K°q y , where 

K a = — ; . This eddy diffusion coefficient is used in the transport equation 

w k 2 ( u-cy+8 

for all tracers. When the zonal wind (u ) equation is solved, he compute the EP flux 
divergence as the eddy potential vorticity directly from the solution of the wave equations 
in order to include the contributions not only from the wave breaking but also from other 
processes (e.g. thermal damping of the wave). Kinnersley and Harwood (1993) and 
Kinnersley (1995, 1996) modified Garcia's method in an isentropic two-dimensional 
model by solving the eddy Ertel potential vorticity equation for the three longest waves. 

A wave-breaking criterion similar to the one used in Garcia (1991) was applied to each of 
the three waves. A simple procedure was used to find the damping rate (8 m ) which is 

large enough to stop the growth of the wave where its breaking criterion is met. The 
resulting eddy diffusion coefficient is the sum of the contribution from 3 waves, 

k k _ f S m W 

” fckJiu-cJ+Sj • 

Our approach differs from the above approach in that, the eddy diffusion 
coefficients for each chemical tracer is derived directory from the time dependent wave 
transport equation for the tracer and are calculated using wave fields simulated in a semi- 
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spectral dynamical model. Therefore, it is not necessary to make the assumption that the 
eddy diffusion coefficients derived based on the flux-gradient relationship of PV can be 
used for chemical tracer. 

3. Model description 

The overall model has three modules, a semi-spectral dynamical model, a 
radiative transfer code, and a 2-D CTM. The model has 9.5° resolution in latitude (19 
boxes from pole to pole), ~3.5 km resolution in altitude (half scale height), and the three 
longest zonal waves are calculated in the dynamical module. The 2-D CTM has 17 
vertical layers (from the surface to ~60 km) while the dynamical module includes seven 
extra layers to place the upper boundary at ~84 km. 

The three-wave model differs from most interactive 2-D models in that it uses a 
highly truncated three-dimensional dynamical module to calculate the zonal mean and 
wave fields. In the current version of the model, only the three longest planetary waves 
are included. The interactive coupling between the dynamics and the distribution of 
atmospheric constituents is achieved when the model-calculated O 3 is used to compute 
the radiative heating rates. The model flow chart is shown in Fig. 1. Interactions 
between different processes and feedbacks can be identified by following arrows in that 
figure. The data are exchanged between the dynamical module, the radiative transfer 
module and the chemistry transport module every ten days. 

3.1 Dynamical Module 
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The dynamical module is adopted from Schneider and Geller (1984). It is 
formulated in the primitive equation system and the vertical coordinate is z = H*ln(po/p). 
The zonal mean equations are 


d t u + tan0j + w9 z u =-E v + Dfj ~P R u, (1) 

9,v + u(f + — tan G) = -dJ> - E v + D v ~P R v , (2) 

a 

9& + v9j z + w9j z + wN 2 =±J-E t + D t , ( 3 ) 

H 

—^—9{vcos8) + — 9 z (p 0 w) = 0. (4) 

COS 8 On 


Here 9 denotes the partial derivative, e. g. d t is the time derivative, u, v and w are the 
components of the velocity, 0 is the geopotential, and <j> z = d z <p , which is proportional 
to the temperature, 6 is the latitude, a is the radius of the Earth, p 0 is the reference air 

density, f is Coriolis parameter, N is the Brunt- Vaisalla frequency, H is the scale height, 
taken to be 7 km, K=2/l, is the Rayleigh friction coefficient which is a function of 
altitude. J is the radiative heating. E u , E v and Et are eddy flux divergence , 


E u = — d (« V cos 8) + — d z (p 0 uw') - tan 8 , (5) 

cos 0 p 0 a 

E v = —^—9 (v'v'cos8) + —d (p 0 v'w') + ^-tan6 , (6) 

cosd p 0 a 

E T =-^—d y (7r z cosd) + — d z (p Q W z ) ■ (7) 

COS 8 Pn 


The equations for deviation from the zonal mean (eddy), denoted by a prime, are 

9/ + u9 x u' + v'9 u — fv' - UV - + - U — tan 8 = -9 x <p' -W u +D u .~ P R u' (8) 

a 

9/ + u9y + v'<9 v + w'9 z v + fu' + 2 tan 8 = -9 J>' -W v + D v .- p R v' (9) 

a 
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( 10 ) 


d,<t>' z + ud x <p' + v'd y (t> z + w'd z <t> z + w'N 2 = -J'-W T + D r 

H 

d x u + —1— d ( v'cos 0) + — d ( p 0 w ') = 0 . (11) 

cos 9 p 0 


where W u , W v and W T are the terms of second order in the eddy equations, which 
represent the wave-wave interactions. 


W u = d x (u ' 2 ) + d ((mV - mV) cos 9) + 


COS0 


+ — d z (p 0 (mV - mV)) - 

P o 


mV -mV 


a 


tan 9 


( 12 ) 


= d x (u'v') + — l — d f(v' 2 -v' 2 )cos0) + 

COS ^ _ 

+— d ,(p 0 (vV — vV)) + — — tan0 , (13) 

Po ^ 


= <9,(m'0;)+^-^((v>;- v>;)cos0+ 
cos 6/ 

+— ^(p 0 (w>;-v^)). d4) 

Po 


D u , D v , D T in Eqns. (l)-(3) and , D v . , D r in Eqns. (8)-(10) represent sub-grid scale 

diffusion and a biharmonic diffusion is used along the latitude direction with time scale 
of 14 days for a 2 Ay-wave and 200 days for a 4 Ay- wave for D u , D v , D r (Schneider 

and Geller, 1984). 


In the stratosphere only a few planetary scale waves prevail according to the wave 
propagating theory (Chamey and Drazin, 1961). Therefore, we keep the three largest 
zonal waves, i.e. wavenumber 1, 2 and 3 for all the variables in Eqns. (8)-(l 1) using 
Fourier expansion around the latitude circle, 
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¥' = ^¥ m (t,y,z)e ik ^ , 

m = \ 


where y/' = (w',v',w',0') , k m = — — — , m is the wave number, 0 is the latitude 

a ■ cos 8 

and a is the radius of the earth. The Eqns. (8)-( 11) then become the wave equations for 
the Fourier coefficients y/ m 's (See Schneider and Geller, 1984 for the details). The y/ m 's 

are complex numbers. Their magnitude measures wave amplitude and their argument is 
related to the wave phase angle. When expanded in Fourier series, the wave-wave 
interactions terms ( W u , W v , W T ) should contain terms in the form of > 

m,n - 1 

where y/ m , £„ can take on u m , v m , w m or <t> zm . Since we are only solving for the first 3 

waves, the summation has to be truncated at m, n = 3. We add to the right hand side of 
the wave form of Eqns. (8)-(10) a linear damping term (- y ss u m ,~ y ss v m , and - y ss <j> Zm 

respectively) to parameterize the wave decay due to the wave- wave interactions with 
smaller unresolved waves (wavenumber m > 3). The same damping rate is used for u , 
v' and 0' . The damping time constant, r„ = l/y„ , is an adjustable parameter. Our 

results show that the value of 60 days gives a reasonable calculated distribution of wave 
amplitudes for geopotential height. 

It should be pointed out that a severely truncated model may result in kinetic 
energy build-up on the shortest wave (O'Brien and Branscome, 1989; Zhou and Stone, 
1993). This does not occur in our three-wave model because the model has relatively 
fine resolution in the meridional direction (9.5°) and large sub-grid scale diffusion. 

The model's vertical resolution in the troposphere is too crude to generate 
stationary waves via surface topography. We specify the lower boundary condition for 
wave amplitudes and phases from a set of climatology data. The data are extracted from 
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an archive of monthly mean grid point analyses of meteorological data (See the appendix 
for details). We only specify the first three wave components of the monthly mean 
geopotential height ( $ ) at 850 mb, the lower boundary of the model. The amplitudes of 
wave geopotential height at the lower boundary are shown in Fig. 2 for January and July 
respectively. The monthly mean data are interpolated to each time step and at the bottom 
layer of the model when the wave equations are integrated. The zonal mean and wave 
equations are solved with a time step of 30 minutes. 


3.2 Radiative Transfer Module 

The radiative heating in the zonal mean thermodynamical Eqn. (3) is calculated 
every 10 days using a narrow -band radiative transfer scheme described by Wang and 
Ryan (1983). In the eddy thermodynamic Eqn. (10), we use a simple Newtonian cooling 
form to represent J' , i.e. J' = —OCT'. The Newtonian cooling coefficient OC(6,z) is 
calculated from the zonal mean radiative heating using 

J(T + AT,d,z)-J(T,0,z) 
a(tf,z) = — . 


with AT taken to be 0.5 °K. 

Radiative heating is not the dominant component of the heating rate in the 
troposphere. The release of latent heat during convection is more important. The model 
does not directly calculate the latent heat release. Instead, the Cunnold (1975) 
parameterization scheme is used to include all heating sources by introducing in the 
energy equation a term that represents relaxation towards a prescribed temperature field 
with a prescribed time constant. The time constant is specified as a function of altitude 



only, which is 10 days at the surface, increased to 20 days at 10 km and kept equal to 20 
day up to the tropopause. This constant is also used as the Newtonian cooling time 
constant in the troposphere to calculate J' . 

3.3 Chemistry Transport Module 


As in most 2-D chemistry transport models, the grouping technique is used in the 
chemistry transport module of the three-wave model. All the fast chemical species are 
assumed in the local photochemical equilibrium state and are not transported 
individually. They are grouped into families. Each family species is transported with the 
appropriate net production and removal rate for the group. Our model includes 19 slow 
species (N 2 0, CH 4 , C 2 H 6 , H 2 , CO, CH 3 C1, CCI4, CFCI3, CF 2 C1 2 , CH3CCI3, CF 2 0, 
CH 3 Br, CHC1F 2 , C 2 C1 3 F 3 , CBrClF 2 , CBrF 3 , O3, HNO3, H 2 0) and five family species 
(NOy, CLy, Br y , Fy, NO z =NOy - HNO3). 


The zonal mean transport equation for a slow species or a family species is 


d,X + v *dyX + w 'd z X = S-yx~VF 


ft 

dz 

(15) 


~d z 

Pa 


-%v'# 

N 2 Vz 


3x 


1 


jdy 


cos 0^ y 


cos 0 
N 2 


*'<K 


where % I s the zonal-mean mixing ratio of chemical species, S denotes the zonal- 
mean production rate, yx is the zonal-mean removal rate. Here we have neglected the 


second-order term y' ■ x' , which is assumed to be smaller than y % ■ 



In Eqn. (15), v* and w* are velocities of the residual mean meridional 
circulation, defined as 


Po 


— * — 1 -i / COS 0 ~TT}\ 

W =W+ X^v(^72“ V ^)> 

cos# N 


( 16 ) 


and F is the eddy flux, 

F = (17) 

We could solve Eqn. (15) using values for F calculated from a numerical solution of the 
eddy equation for %' : 

d,X' + udX + v'dyX + w'djc = -7x'-r'X- YssX' + S'. (18) 

Note that we added a small scale mixing term - y ss x' in the equation. This term 
represents damping of zonal asymmetry in the tracer distribution by short-waves working 
in connection with planetary scale waves (cf. Matsuno (1980)). Eqn. (18) can, in 
principle, be solved in spectral form to obtain eddy tracer distribution X - However, in 
the dynamics module we don't directly calculate w' , which is needed in order to solve 
Eqn. (18). Also there are several numerical difficulties associated with this approach. 
Apart from the well known problems related to the spectral solutions of the continuity 
equation (Allen et ah, 1991), including y' ■ x and S' terms will couple the eddy 
transport equations for different species. 



In this study, we take the alternative approach by recognizing that if we ignore the 
y' ■ x and S' terms in Eqn. (18), Eqn. (17) can be written in the form (See e.g. Strobel’s 
derivation using a linear wave (1981)), 

d,X + v‘d y x + wd z X = S-fx + V(K- V£), (19) 

where K , the diffusion tensor, is approximately symmetric (See Tung 1984; Tung et al. 
1985). Tung (1984) argued that the eddy transport in the stratosphere due to planetary 
waves is dominated by the mixing along isentropic surfaces. Because the isentropic 
surfaces are nearly horizontal, it is usually assumed that K yy is equal to the diffusion 
coefficient along the isentropic surface and Ky Z , K Z y can be calculated as the projections 
of K yy from the isentropic surface back to the pressure coordinates. 

In order to solve Eqn. (19), the chemistry scheme from the AER two-dimensional 
chemical transport model (Ko et al. 1984, 1985, 1989; Weisenstein et al. 1993) is used to 
calculate the chemical source term 5 and loss rates / . The same numerical scheme 
described in the 2-D CTM is used to solve the transport equation for trace gases with a 
time step of 8 hours. The residual circulation ( v * , vv* ) is calculated according Eqns. 

(16) and (17), and the v'<p' z term required is computed at each (dynamical) time step of 
half an hour and averaged over 10 days in the dynamical module. At the end of each 10 
day period, they are sent, along with a set of v and vv , to the chemistry transport 
module. The values of v* and vv* are then formed using Eqns. (15) and (16). The 
values used in each eight hour time step are obtained by interpolating between v* and 
vv* calculated 10 days apart. 

The formulas for calculating the stratospheric Kyy, which describes mixing along 
isentropic surfaces, from dynamical wave components, chemical loss rates and small 



scale mixing are given in the next section. A constant K zz of 10 rn^/s in the troposphere, 
0.1 m^/s in the stratosphere, and 1 m^/s in the mesosphere is used to parameterize the 
rapid mixing by convections and turbulence in the troposphere and gravity wave breaking 
in the mesosphere. To represent the fast turbulent mixing in the troposphere, a constant 
value of 1.5 x 10 6 m 2 /s is used for K yy in the troposphere except along the tropopause, 
where the values for K yy are adjusted to represent the exchange between the stratosphere 
and the troposphere (Shia et al., 1993). These diffusion coefficients, along with 
calculated residual mean circulation, are then used in the zonal mean chemical transport 
model to calculate the distribution of the trace gases. 


4. Eddy Diffusion Coefficients K yy 


Assuming that j^'and the perturbations in the velocity, v' and w' can be 

represented by linear waves, with the same longitudinal and time dependencies, i.e. 
exp(-cff + i(kx - cot)). Strobel (1981) showed that 



-VdyX ~ 

( y -a )- i{co - ku) ’ 


( 20 ) 


and K yy is the coefficient of d v in v'%' , given as 



y -a 


(y-a)2 +(co-ku)2 


vV. 


( 21 ) 


We can generalize the expression for K yy for the multi-wave case with small scale 
mixing if we assume that %' , V and w' share the same x and t dependencies, 

f'(t,x,y,z) = Re(£f m (t,y,z)e ,kmX ) = Re(^/ m (y,z)e' ,(<0 "' , “" ) 'e ,i " j: ), 
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m=] 


( 22 ) 



where CC m (wave decaying rate) , and (D m (wave frequency) are functions of latitude and 
altitude only, is negative for a growing wave and is positive for a decaying wave. 

The expression for %' analogous to Eqn. (20) in the multi- wave case and with 

small scale mixing is given by 


~ v J y X ~ wAX 


(23) 


2 N 

By examining v'x' = — Rc(^^ m (r,y,z)v* (t,y,z)), we obtain an expression for K yy as 


m = 1 


following, 


K ,o, _ y , (r+r^-^J-Mw)! 2 , 

* sV.-^+ca.-7-rJ 2 

For diagnostic purpose, we define 


(24) 


TStot zrtr ■ T/'ch , tsss 

yy ~~ yy yy yy* 9 


where 




~«m •|V W (^) ; .Z)| 2 


i (G) m -k m uy+(a m -y ch -Yj 




K = I(; 


yc/,k( f »y»^r 


( 0) m -k m u) 2 +(a m -Y ch -Yss )' 


m = 1 




(25) 


*£ = X< 


yss-\ v Jyy ^)\ 2 


^ (a m -k m u) +(a m -Y ch -Yj 


-)• 


The truncating at wave number 3 is justified by assuming that the direct 
contribution from the smaller unresolved (m>3) waves to the global transport is 



negligible compared with that from the planetary waves. However, the motions 
associated with wave number m > 3 create small scale mixing which can be coupled 
with planetary waves to generate global scale mixing. The y a is the damping rate for 

small scale mixing. We assume that this small scale mixing is related to the small scale 
mixing in Eqns. (8)-( 10) in the dynamics model and the same time constant, r ss = 1 !y a , 

of 60 days also should be used here for tracer mixing. This number is quite close to 
Matsuno's (1980) estimated value for small scale mixing of 10 7 sec. 

The K‘ y r y defined in Eqn. (25) is the transient eddy coefficients. It is 
approximately proportional to the wave growing rate, -a m , which is species 

independent, and represents the mixing of inert gas by transient waves. This is consistent 
with the conclusion in Plumb (1979) and Pyle and Rogers (1980), that steady waves ( a m 
= 0) do not transport chemically inactive tracer. A growing wave ( a m < 0) implies 
positive K" y , giving rise to down-gradient transport. This kind of transport during wave 

growing periods was confirmed by observations (Leovy et al., 1985). After the wave 
reaches its peak, it decays ( a m > 0). For decaying period, K' )y is less than zero and 

produces a up-gradient transport. The up-gradient transport in decaying period 
compensates the down-gradient transport in the precedent decaying period and generates 
little net transport. However, net transport results when the local tracer concentrations 

are changed by chemical reactions or by small scale mixing. They are represented by 
K* and AT" , and both are always down-gradient (diffusive). 


In order to determine K yy , we need the wave frequency, a> m and wave decaying 
rate a m . They can be determined by 


d,v„(f.y,z) _ 

v m (t,y,z) 


-ICO. 
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using v m calculated in the dynamics model. In doing so, we implicitly assume that one 
planetary wave is dominant in the zonal asymmetry with wave number equal to m . This 
is probably true in the region where the planetary waves are strong and generate large 

eddy mixing, which is what we are parameterizing. In the stratosphere, the typical time 

scale is about 10 days for ft) and several weeks for a m . 

3 m m 

The diffusion coefficients, , K* , and K™ are calculated in the dynamical 

module for every 30 minutes time step using Eqns. (25) and (26). The mean values of the 
eddy diffusion coefficients over the 10 day period are used in the transport equation for 
trace gases. 


5. Model results 

5.1 Dynamical Variables 

The calculated global distributions for mean zonal wind and zonally averaged 
temperature for January and July are shown in Fig. 3. The magnitudes and locations of 
the polar jets in the winter hemisphere are in good agreement with observations (e. g., 
Barnett and Comey, 1985). The seasonal changes of subtropical tropospheric jets are 
also well simulated. However, the upper easterly in the summer hemisphere are weaker 
than observations. The winter zonal mean temperature in the polar lower stratosphere is 
more than 10° colder than observations. This problem is common in models with 
moderate resolution. 

The latitude-height section of wave amplitude of wave 1, 2 and 3 calculated by 
the model and from the Climate Diagnostics Data Base (CDDB; see the Appendix) data 
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for January and July are plotted in Fig. 4. In general, the model simulates the wave 
amplitudes well except that the amplitude of wave 1 is much smaller than the 
observations during winter in Northern Hemisphere. The parameterization of the heating 
rates in the troposphere limits the interaction between the troposphere and the 
stratosphere and may distort the wave propagation. Since the stratospheric waves 
propagate from the troposphere, the very strong hemispheric asymmetry of eddy activity 
in the stratosphere is mainly attributed to the specified lower boundary conditions. For 
the CDDB data, which we have used to prescribe geopotential height, waves are much 
weaker in the southern hemisphere winter than in the northern hemisphere winter, 
especially for wavenumbers 2 and 3 (Fig. 2). 

5.2 Residual Mean Circulation 

We calculate the residual mean circulation from the Eulerian mean circulation and 
eddy heat fluxes explicitly, according to Eqns. (16) and (17). The stream functions of 
residual mean meridional circulation are shown in Fig. 5 for four seasons. The stream 
patterns are very similar to those used in the TEM models or from the diabatic circulation 
(e.g., Hitchman and Brasseur, 1988; Ko et al., 1985). The upwelling motion occurs in the 
tropics and summer subtropics, and downward motion in the winter hemisphere. For 
spring and fall, the downward motion occurs in both hemispheres. These are the features 
of diabatic circulation. There has been some concern about the accuracy of the residual 
circulation calculated directly from Eqns. (16) and (17) due to the cancellation of the two 
terms on the r.h.s. of the equations. It is apparent from our model results that this is not a 
serious problem to affect the accuracy. 


5.3 Eddy Diffusion Coefficients 
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The calculated total eddy diffusion coefficient, K ™ for ozone and eddy diffusion 
coefficient for an inert tracer, which is K yy + K”, are plotted in Figs. 6 and 7 respectively. 

Their difference is the chemical eddy due to the local chemical reactions which ozone 
participates. Both coefficients are smaller than 4 x 10^ m^/sec in most areas and all 

seasons except in the northern middle latitudes during winter. The tropical stratosphere 
has very small eddy diffusion year around. Negative values of K ™ are calculated for 

ozone only during spring in northern upper midlatitudes. The magnitude of the negative 
values of K™ is quite small. The time step in the numerical scheme for solving Eqn. 

(20) is made smaller when necessary to avoid negative mixing ratio in the presence of the 
negative K yy . As expected, K™ for ozone is larger than K™ for an inert tracer . Also, 

K™ for long-lived species, e. g. methane (not shown), is very close to K™ for inert 
tracer. 


As a comparison, the values for K yy , , which is derived from the flux-gradient 

relationship of the potential vorticity, are calculated using the model calculated winds. 
The quasi-geostrophic potential vorticity and its eddy are calculated using the following 
expressions. 

q = f + d x v- d y u +^d z (£j&) (30) 

Po 

q' = d x v' - d ll' (31) 

y Po 


Here / is the Coriolis parameter, p 0 is the reference air density, which is the function 

of z only, and 0 is the potential temperature. Using the least-square method 
(Newman et al., 1988), the values of V ■ q' and q are used to derive K ^ . 



( 32 ) 



y' • q' • q, 
[q y f 


The direct application of Eqn. (32) will give large negative values of in some 
regions where q is close to zero and the eddy flux is up-gradient. Therefore, a lower 
limit of zero is applied when K ^ is used in the transport calculations. 


Model calculated K yy shown in Fig. 8 shares the same features with that in 

Newman et al. (1988), such as north-south asymmetry, the seasonal variation. But the 
details in two set of K q yy are quite different due to the different wind fields used to 

calculate the potential vorticity. We are more interested in the differences between the 
model calculated and K™ for ozone and other tracers, because both are calculated 

using the quantities generated in the same model. As displayed in Figs. 6 , 7 and 8 , K yy 
is larger than the total eddy diffusion coefficient, K'° y for ozone and for inert tracer. 

One feature is shared by all these eddy diffusion coefficients: small diffusion in 
subtropical stratosphere. 

The sign of the transient eddy diffusion coefficient, K \ * , changes frequently, and 
K yy is always positive. Their magnitudes are comparable so the total eddy diffusion 

coefficients are mostly positive (Fig. 6 and 7). 


5.4 Long-lived Species 

N 2 O, CH 4 are good tracers because their chemical reactions are slow in the 
stratosphere and their sources are at the surface. Their distributions in the stratosphere 
are indicative of transport processes. To evaluate how realistic the calculated residual 
mean circulation and eddy diffusion are, we compare the model results of N 2 O and CH 4 
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with the UARS/CLAES (Kummer et al., 1993) and UARS/HALOE (Tuck et al., 1993) 
data. The data show that the mixing ratios of N2O and CH4 have two important features. 

The first is a "double peak" structure in the upper stratosphere during the period from 
February to May. This feature is attributed to a strong tropical meridional circulation 
(Choi and Holton, 1991). The second is that the tropical maxima tend to lean toward the 
summer hemisphere. In the upper stratosphere the maximum location is about 30° away 
from the equator in solstices. From December to May the ridge leans toward south, while 
from June to November it leans toward to north. The model calculated mixing ratios of 
N2O and CH4 are shown in Figs. 9 and 10 respectively. The "double peak" structure 

only vaguely shows up in the model simulations in the upper stratosphere during the 
northern winter. However, the "leaning ridge" feature is apparent in the summer 
hemisphere for both species and the simulated positions of the "leaning ridge" are close 
to where they are observed. 

Simulated distributions for both CH4 and N2O have steep "slopes" and large 
meridional gradients in the extratropics. This feature is also seen in recent observations 
of N2O and CH4 from UARS/CLAES (Fig. 11). It is noted that these features have not 

been successfully simulated by some two-dimensional models, probably due to large 
eddy diffusion used in these models. As a demonstration, the simulated CH4 
distributions for four seasons using K* is plotted in Fig. 12. They are more 'flat' 

compared with Fig. 11. 

Because both N2O and CH4 are long-lived constituents, they share the same 
mixing surfaces in the lower stratosphere. According to Plumb and Ko (1992), the 
correlation between two long-lived species is compact if the horizontal diffusion 
dominates the transport process. The correlation between N2O and CH4 in the northern 

winter tropics and extratropics is shown in Fig. 13. Two compact groups of points are 
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displayed there, one, which is approximately linear, in high latitudes (poleward of 57° N) 
and one in tropics and subtropics. The same feature also can be found in the Southern 
Hemisphere and in other seasons. This can be interpreted as the existence of the tropica] 
pipe (Plumb, 1996) in our model. 

The lifetime calculated from the model results is 109 yr for N 2 O, 57 yr for F-l 1 
and 95 yr for F-l 2. Compared with the values of 132 yr for N 2 O, 57 yr for F-l 1 and 1 12 

for F-l 2 from AER 2-D CTM without the tropical pipe ( For the detail discussion of the 
tropical pipe, see Plumb, 1996); and of 103 yr for N 20 , 45 yr for F-l 1 and 86 for F-l 2 

with the tropical pipe (Weisenstein, private communication), the three-wave model 
results are between that of these two CTMs. Because all three models use the same 
photochemistry, the differences are due to the different residual mean circulations and 
eddy diffusion coefficients. 

5.5 Ozone 

The O 3 column abundances from the model results and from TOMS data are 
plotted in Fig. 14 as a function of month and latitude. The ozone column abundance from 
the baseline model results (Fig. 14 a) shows a strong hemispheric asymmetry. The 
sources of the asymmetry are the differences of the lower boundary condition for wave 
generation and the Cunnold (1975) parameterization of tropospheric heating rate in the 
two hemispheres. These differences in the model inputs generate asymmetry in the 
model calculated eddy diffusion, the mean meridional transport circulation and the tracer 
distribution. The general features of the O 3 column distribution are similar to the 

observations from Total Ozone Mapping Spectrometer (TOMS), which are shown in Fig. 
14 b. The modeled spring maximum in both hemispheres appears in the right places with 
roughly right amplitudes but a little later in time. The calculated ozone column 



abundance in the northern mid- and high-latitudes during late summer and fall is about 
40-60 Dobson Unit (D. U.) higher than the observations. 

It is interesting that the ozone distribution, (Fig. 14 c), calculated using the eddy 
diffusion coefficient of inert gas ( K™ with y ch =0) for all tracers is quite similar to that 

using the total eddy diffusion coefficient, (Fig. 14 a). This means that the contribution 

from the chemical eddy diffusion is not essential in our model. In contrast, the ozone 
distribution calculated using , (Fig. 14 d), is quite different. Due to larger values of 

, the ozone maximum in the northern polar spring is 40 D. U. less than observations, 

but the ozone column abundance in the same region during fall is improved, compared 
with baseline simulation with K . 

5.6 Correlation Between Ozone and NOy 

Observations have shown that mixing ratios of O 3 and NOy are well correlated in 
the lower stratosphere (Murphy et al., 1993; Fahey et al., 1996). The ratio of NOy to O 3 
is nearly independent of altitude from the tropopause to above 20 km. The ratio shows 
strong meridional gradients in the subtropics. Previous two-dimensional models (e.g., Ko 
et al., 1993; Solomon and Garcia, 1984) have simulated general features of the vertical 
profile of NOy /03 but not the gradient in the lower stratosphere between tropics and mid- 
latitudes. The mixing ratio of NOy, O 3 and their ratio calculated by this model are 
shown in Fig. 15. The feature of NOy /03 is very similar to the observation, especially, 
the strong subtropical meridional gradients are reproduced. 


6. Discussions 
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The model reproduces many observed characteristics of dynamics and tracer 
distributions. The waves are forced by the assigned boundary condition at 850 mb and 
propagate through the troposphere into the stratosphere. The simulated zonal mean fields 
and eddies give rise to a reasonable set of residual mean circulation, temperature and 
eddy diffusion coefficients to the chemistry transport module, which generates 
distributions for ozone and other long-lived species with desired features. The 
distributions of long-lived species like N2O and CH4 demonstrate a subtropical transport 

barrier in the lower stratosphere due to the upward residual circulation in tropics and the 
small eddy diffusion coefficient in the subtropics. The calculated values for K™ are 

smaller than 4 x 10^ (m^/sec) from 30° N to 30° S in the lower stratosphere all year 
around (Fig. 6). Therefore, the tropics is isolated with little exchange between tropics 
and middle latitudes and acts like a pipe (Plumb, 1996). This also limits the direct 
exchange between the two hemispheres in the lower stratosphere. The calculated 
correlation between N2O and CH4 does distinguish the tropical area from the high 

latitudes. 

The total eddy diffusion coefficients for inert tracers ( K™ with f ch =0, shown in 

Fig. 7), are smaller than the eddy diffusion coefficients derived based on the flux-gradient 
relationship of potential vorticity, K^, in our model (Fig. 8). We think that a part of the 

difference is due to the dissipation of the potential vorticity in our model. In the current 
version of the model, the dissipation consists of two parts, radiative cooling and friction. 
With the dissipation, the potential vorticity behaves like a chemically active tracer. 
Therefore the eddy diffusion coefficients derived from the PV flux-gradient relation 
include both transient and "chemical" eddy contributions for the PV. Our model 
calculations show that the dissipation does have a significant contribution to the eddy 
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diffusion coefficients for PV. Otherwise K q y would be close to the total eddy diffusion 
coefficients for inert tracers. 

Most 2-D CTMs neglect the chemical eddy diffusion and use the same eddy 
diffusion coefficients for all trace gases. Our model results demonstrate that the chemical 
eddy diffusion is unimportant in simulating the distribution for most slow species. The 
mixing ratio of most slow species including ozone, N 2 O, CH 4 , C 2 H 6 , H 2 , CO, CF 2 O, 
NO y and Cl y change only less than 10% in most regions and all seasons when chemical 
eddy diffusion, K‘ y , is ignored. We think this is the result of competition between the 

transport and photochemistry. In the region where the chemical eddy diffusion 
coefficient, which is proportional to the local chemical loss rate, f , is large, the local 
lifetime ( 1 / y ) of this species must be short and the distribution is mainly determined by 
the photochemistry not by the transport. In the region where the chemical eddy diffusion 
coefficient is small the eddy mixing is dominated by the transient eddy and small scale 
mixing. Therefore, in either case the chemical eddy transport is not essential. 

Despite the general agreement between the model simulation and the 
observations, several aspects of the model results need to be improved. With finer 
vertical resolution we expect better wave simulation and meridional energy transport to 
make the lower stratosphere poles warmer during winter. Better representation of the 
tropospheric heating rates will improve the tropospheric dynamics, the interactions 
between the stratosphere and troposphere, and the simulation of wave propagation. The 
calculated ozone column abundance (Fig. 14 a) in the northern high latitudes during fall 
is about 40 to 60 Dobson Unit higher than observation (Fig. 14 b). We expect that the 
future improvements of the model (e. g. increasing the number of waves resolved in the 
model; using finer vertical resolution) will alleviate this problem. 
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Although the model is built on the primitive equations and other well established 
equations in physics and chemistry, in order to make the model manageable certain 
approximations and assumptions have to be applied. The followings are some of the 
significant assumptions to the model formalism: 

(1) the truncation at m = 3 in the wave expansion and adding the linear damping term in 
the wave equations to represent dissipation caused by the wave-wave interactions with 
small waves (m > 3); 

(2) parameterizations for the tropospheric heating rates, gravity wave dragging (Rayleigh 
friction), sub-grid diffusion in the latitude direction; 

(3) using family grouping and assuming y' • %’ « 7 ‘ X solving zonal mean tracer 
transport equations; 

(4) ignoring S' and y' • % in the eddy transport equation for tracer; 

(5) assuming that the tracer waves, and the velocity waves, V m 

have the same decaying rates in calculating K ; 

(6) using a linear damping term in the tracer transport equation to represent mixing by 
small scale waves. 

Most of the approximations and assumptions we used are generally accepted in 
literature or reasonable. Some of them can be verified in principle as satellite data or 3-D 
model results become available. 
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APPENDIX 


The data used for the model lower boundary condition and for diagnosis were 
extracted from an archive of monthly mean grid point analyses of meteorological data, 
known as the Climate Diagnostics Data Base (CDDB) (Rasmussen and Arkin, 1985), 
based on the U. S. National Meteorological Center global analyses, which are generated 
at 0000 GMT and 1200 GMT daily. The original data were on a 2.5° by 2.5° grid, at nine 
standard levels in the atmosphere (1000 mb, 850 mb, 700 mb, 500 mb, 300 mb, 250 mb, 
200 mb, 100 mb and 50 mb). The data for use in the model consisted of the first seven 
complex Fourier components interpolated onto the eighteen edges of the 19 bands of 
latitude used by the model, over the same nine levels. 

The original data set spanned the time period from October 1978 to March 1990, 
and contained values for the zonal component of the wind velocity (u), the meridional 
component (v), temperature (T), geopotential height ( <f > ), vertical velocity, specific 
humidity, vorticity and several cross-products. The first step in the extraction 
process was to create a subset of data from October 1978 to September 1989, which 
contained u, v, T and <j > , at each of the of the nine pressure levels. The climatological 
considerations of some of these variables taken from the CDDB for a similar period, the 
1980-89 decade, including estimates of the seasonal and interannual signals, are given by 
Salstein (1995). 

After the spatial series were created, Fourier transforms in space were calculated 
at each latitude, for each month, variable, and level. The first seven complex harmonics 
of each transform generated were stored. A climatology was created for the October 
1978 to September 1989 period and is used in this study. 
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Finally, because the model used nineteen bands of latitude, whose edges were 
defined by eighteen equally spaced latitude circles, and the model required input values 
along these edges, area weighted averages of the complex Fourier transforms were 
calculated around each of these latitude circles. 
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Figure captions 


Fig. 1 The flow chart of the three-wave interactive model. The arrows indicate the 
direction of data flow. The loops formed by arrows represent feedbacks. 

Fig. 2 The prescribed wave amplitude in the model's lower boundary, for (a). January and 
(b). July respectively. The solid lines is for wave 1, dash line for wave 2, and dot line for 
wave 3. 

Fig. 3 The model calculated mean zonal wind, u , and temperature, T , (a), u in 
January, (b). T in January, (c). u in July, and (d). T in July. The contour interval is 
10 m/s for u and 10° for T. 

Fig. 4 Comparison of the model results and the observations for the wave amplitude of 
geopotential height, (a) January; (b) July. The contour increments are 40 m, 20 m, and 10 
m for wave 1, 2, and 3 respectively. 

Fig. 5 The calculated streamfunction of the residual circulation, y / , for (a) January, (b) 
April, (c) July, and (d) October. The contours plotted are -2,-1, -0.5, -0.2, -0.1, -0.05, - 
0.02, -0.01, -0.005, 0, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2 x (10 8 km 2 /s). The 
values correspond the meridional mass flux scaled by (A p s ), where A is the radius of 
the Earth and p s is the surface air density. 

Fig. 6 The total eddy diffusion coefficient, K ™ , (10^ m^/sec) for ozone, in (a) January, 

(b) April, (c) July, and (d) October. The contour levels are 0, .05, 0.1, 0.2, 0.4 , 0.8, 1.0 
and 1.5. The areas with negative values are shaded. The calculated K‘°‘ is plotted for 
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whole model region. In the model, a constant value of K yy is added in the troposphere to 
simulate rapid mixing in the troposphere. 

Fig. 7 Same as Fig. 6 , but for inert tracer ( K™ with y ch =0). 

Fig. 8 Same as Fig. 6 but for eddy diffusion coefficient derived based on PV flux- 
gradient relation, K yy . The comparisons of K yy with the total eddy diffusion 

coefficient, K™ , for ozone shown in Fig. 6 or for inert tracer shown in Fig. 7 should be 
limited in the stratosphere. 

Fig. 9 Latitude-altitude cross section of calculated N 2 O mixing ratio (ppbv) for (a) 
January, (b) April, (c) July, and (d) October. The contour levels are 5, 10, 25, 50, 100, 
150, 200,250, and 300 ppbv. 

Fig. 10 Latitude-altitude cross section of calculated CH 4 mixing ratio (ppmv) for (a) 
January, (b) April, (c) July, and (d) October. The contour increment is 0.2 ppmv. 

Fig. 1 1 The Monthly averaged zonal mean of (a) UARS/CLAES N 2 O for March of 1993, 
and (b) UARS/CLAES CH 4 for March of 1993. Archived data are from UARS CD- 
ROM Vol. 7. 

Fig. 12 Same as Fig. 10, but using eddy diffusion coefficient derived based on PV flux- 
gradient relation, K“L , instead of K yy . 

Fig. 13 The correlation of N 2 O and CH 4 mixing ratio (a) in the northern mid-high 
latitudes, and (b) in the tropics. A straight line is 
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added in each panel to show that the transport characters in the tropics and in the high 
latitudes are different (Plumb and Ko, 1992). 


Fig. 14 Latitude-time cross section of ozone column abundance (D. U.) from (a), model 

calculation, (b). TOMS measurements, (c). same as (a) but using the eddy diffusion 
coefficient for inert tracer ( AT"' with f ch =0), for all trace gases, and (d). same as (a) but 

using eddy diffusion coefficient derived based on PV flux-gradient relation, , for all 
trace gases instead of AT"' . The contour increment is 20 D. U. 

Fig. 15 The mixing ratio of (a) NOy and (b) O3, as well as (c) their ratio, NOy/03 along 
with the AAOE, STEP, and AASE data (Murphy et al., 1993) in the lower stratosphere 
(64 mb), in January. 
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Fig. 1 The flow chart of the three-wave interactive model. The 
arrows indicate the direction of data flow. The loops formed by 

drrni \rc fW>rlhprV<; 






Wave Amplitude (m) Wave Amplitude (m) 


(a) 



(b) 



Latitude 


Fig. 2 The prescribed wave amplitude in the model's lower boundary, 
for (a). January and (b). July respectively. The solid lines is for wave 
1, dash line for wave 2, and dot line for wave 3. 
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Fig. 3 The model calculated mean zonal wind/r, and temperature,? 7 
(a), u in January, (b). f in January, (c). u in July, and (d). f in 
July. The contour interval is 10 m/s for u and 10° for T . 
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Fig. 4 (a). Comparison of the model results and the observations for 
the wave amplitude of geopotential height for January. The contour 
increments are 40 m, 20 m, and 10 m for wave 1, 2, and 3 
respectively. 
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Fig. 4 (b). The same as Fig. 4 (a) but for July. 
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Fig. 5 The calculated streamfunction of the residual circulation,^, 
for (a) January, (b) April, (c) July, and (d) October. The contours 
plotted are -2, -1, -0.5, -0.2, -0.1, -0.05, -0.02, -0.01, -0.005, 0, 0.005, 
0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2x (10-8 kmVs). The values 
correspond the meridional mass flux scaled by (A p,), where A is the 
radius of the Earth and p, is the surface air density. 
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Fig. 6 The total eddy diffusion coefficient, K , (10^ m^/ sec) for 
ozone, in (a) January, (b) April, (c) July, and (d) October. The contour 
levels are 0, .05, 0.1, 0.2, 0.4 , 0.8, 1.0 and 1.5. The areas with 
negative values are shaded. The calculated K™ is plotted for whole 
model region. In the model, a constant value of K„ is added in the 
troposphere to simulate rapid mixing in the troposphere. 
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Fig. 7 Same as Fig. 6, but for inert tracer (K% with y cA =0). 
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Fig. 8 Same as Fig. 6 but for eddy diffusion coefficient derived based 
on PV flux-gradient relation, K* n . The comparisons of with the 

total eddy diffusion coefficient, K" , for ozone shown in Fig. 6 or for 
inert tracer shown in Fig. 7 should be limited in the stratosphere. 
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Fig. 9 Latitude-altitude cross section of calculated N 2 O mixing ratio 
(ppbv) for (a) January, (b) April, (c) July, and (d) October. The 
contour levels are 5, 10, 25, 50, 100, 150, 200,250, and 300 ppbv. 
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Fig. 10 Latitude-altitude cross section of calculated CH4 mixing ratio 
(ppmv) for (a) January, (b) April, (c) July, and (d) October. The 
contour increment is 0.2 ppmv. 
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Fig. 12 Same as Fig. 11, but using eddy diffusion coefficient derived 
based on PV flux-gradient relation, K'„ , instead of K “‘ . 
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Fig. 13 The correlation of N 2 O and CH 4 mixing ratio (a) in the 
northern mid-high latitudes, and (b) in the tropics. A straight line is 
added in each panel to show that the transport characters in the 
tropics and in the high latitudes are different (Plumb and Ko, 1992). 
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Fig. 14 Latitude-time cross section of ozone column abundance (D. U.) 
from (a), model calculation, (b). TOMS measurements, (c). same as (a) 
but using the eddy diffusion coefficient for inert tracer (K% with y ck 
=0), for all trace gases, and (d). same as (a) but using eddy diffusion 
coefficient derived based on PV flux-gradient relation, K'„ , for all 
trace gases instead of K~ . The contour increment is 20 D. U. 
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Fig. 15 The mixing ratio of (a) NOy and (b) O 3 , as well as (c) their 
ratio, NOy/03 along with the AAOE, STEP, and AASE data (Murphy et 
al., 1993) in the lower stratosphere (64 mb), in January. 



